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Abstract. In the ALICE experiment, at forward rapidity (2.5 < y < 4), the 
production of heavy quarkonium states is measured via their fi + fi~ decay channels. 
We present the first measurement of inclusive J/ip production, down to px = 0, from 
Pb-Pb data collected at the LHC at ^/snn == 2.76 TeV. Preliminary results on the 
nuclear modification factor (Raa) and the central to peripheral nuclear modification 
factor (i?cp) show 3/tp suppression with no significant centrality dependence and an 



J* I integrated i?A A 8 ° % = 0.49 ± 0.03(stat.) ± 0.11(syst.) 



Heavy quarkonium states have long been proposed as a sensitive probe of the 



strongly-interacting deconfined medium expected to be formed in the early stages of 
! high-energy heavy-ion collisions [lj. In particular, a measurement of the in- medium 

CO 



dissociation probability of different quarkonium states is expected to provide an 
estimate of the initial temperature of the medium. Suppression of J/ip production 
beyond that expected from cold nuclear matter (CNM) effects (nuclear absorption, 
( ant i-) shadowing, ...) has indeed been observed at SPS and RHIC energies [2J|3llI], but 
several questions are left open. At the LHC, larger suppression might be expected, due 
to the higher initial temperature. However, according to regeneration scenarios |5J, with 
~ 10 times more cc pairs produced in central Pb-Pb collisions compared to RHIC, an 
enhancement of the J/ip production could also be observed. In the ALICE experiment, 
at forward rapidity (2.5 < y < 4), the inclusive production of heavy quarkonium states 
is measured down to pr = via their fi + n~ decay channels in the muon spectrometer, 
as describe in [6]. 

In fall 2010, the LHC delivered the first Pb-Pb collisions at a center of mass energy 
v /snn = 2.76 TeV. ALICE collected data with a minimum bias (MB) trigger, defined as 
the logical AND between signals from the pixel detector (SPD) covering the range \r)\ < 2 
and the two scintillator arrays of the VZERO detector covering the ranges 2.8 < rj < 5.1 
and — 3.7 < r] < — 1.7, in coincidence with two beam pick-up counters, one on each side 
of the interaction region. The centrality of the collision has been determined from the 
amplitude of the VZERO signal [7J . The total data sample available for physics analysis 
amounts to 17 ■ 10 6 MB events. Additional cuts have then been applied to improve the 
purity of the muon sample, in particular by requiring both particles in the pair to be 
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detected in the muon trigger stations. Furthermore, we perform the cuts —4 < rj < —2.5 
and 17.6 < i? a b s < 89 cm, where i? a b s is the radial coordinate of the track at the end 
of the hadronic absorber located in front of the spectrometer, to select muons in the 
geometrical acceptance of the detector. 
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Figure 1. Invariant mass distribution for opposite-sign muon pairs in the centrality 
class 0-10% before (left) and after (right) mixed-event combinatorial-background 
subtraction, with the result of the corresponding fits. 

After these selections, the data sample has been divided into four centrality classes: 
0-10%, 10-20%, 20-40% and 40-80% of the inelastic Pb-Pb cross section. In each of 
these sub-samples, the number ofJ/ip particles has been extracted by using two different 
methods. We have first fitted the opposite-sign dimuon invariant mass distribution 
with a Crystal Ball (CB) function (a gaussian with polynomial tails) to reproduce 
the J/"0 shape and a sum of two exponentials to describe the underlying continuum 
(Fig. [T]-left). Alternatively, we have subtracted the combinatorial background using the 
event-mixing technique and fitted the resulting mass distribution with a CB function 
and an exponential or a straight line to describe the remaining background (Fig. [T]- 
right). The width of the J/ip mass peak depends on the resolution of the spectrometer 
which could, a priori, be affected by the detector occupancy which increases with 
centrality. This effect has been evaluated by embedding simulated J/ip decays to muons 
into real events and no modification has been observed. Such a conclusion has been 
confirmed by a direct measurement of the tracking chamber resolution versus centrality 
using reconstructed tracks. Therefore, the same CB line shape can be used for all 
centrality classes. Several parameters have been tested for the CB tails, either fixing 
them to the values obtained in p-p collisions (where the signal over background ratio 
is more favorable) or in simulations. For each of these choices, the mean and width 
of the gaussian part have been obtained by fitting the integrated mass distribution in 
the centrality range 0-80% and then varying the width by ± 1 standard deviation to 
account for uncertainties (varying the mean has turned out to have negligible effect 
in comparison). The number of 3 ftp particles in each centrality class i (Nj,^), as 
well as the ratios Nj^/N^ 80% , have been determined as the average of the results 
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obtained with the two fitting approaches and the various CB parameterizations, while 
the corresponding systematic uncertainties have been defined as two times the RMS of 
these results, which also approximately corresponds to the maximum difference with 
respect to the mean value. The largest uncertainties have been obtained for the most 
central class and amount to 19% and 12% for N^ 0% and NjJ^ 0% /N^ 80% respectively 

In order to extract the inclusive J /if) yield, N % 3 ,^ has been normalized to the 
number of MB events in the corresponding centrality class and further corrected for 
the branching ratio of the dimuon decay channel and the acceptance times efficiency 
(A x e) of the detector. This latter quantity has been determined from MC simulations. 
The generated J /if) pt and y distributions have been interpolated from existing 
measurements [Sj, including shadowing effects from EKS98 calculations [9J. The 
efficiencies of the muon trigger chambers have been measured directly from data then 
applied to the simulations. For the tracking apparatus, the time-dependent status of 
each electronic channel and their run-by-run evolution have been taken into account 
as well as the residual misalignment of the detection elements. We thus obtained 
a run-averaged value of A x e = 19.4%, with a 7% relative systematic uncertainty. 
The dependence of the tracking efficiency with the detector occupancy has also been 
evaluated using the embedding technique. A small decrease of this efficiency when 
increasing the centrality (—2% in the most central class) has been observed. This 
variation has been confirmed by a measurement of the tracking efficiency performed 
directly from data, and is included in the systematic uncertainties. 

To measure the nuclear modification factors (R AA ), the J /if) yield in the centrality 
class % has been normalized to the inclusive J /if) cross-section measured in p-p collisions 
in the same rapidity domain at the same energy (crj°^ uslve = 3.46 ± 0.13(stat.) ± 
0.32(syst.) ± 0.28(syst.lumi.)/ib) [TO] and scaled by the corresponding nuclear overlap 
function (T AA ) calculated using the Glauber model, while the ratios Nj^/N^ 80% 
have been normalized to the ratios T AA /T AA ~ 80% to extract the central to peripheral 
nuclear modification factors (Rq P ). The systematic uncertainties on the Glauber model 
calculations are 4% and 6% for T^ A 10% and T^ A 10% /T^ 80% respectively. 

The inclusive J /if) Raa are shown in Fig. [2]-left as a function of the average number 
of nucleons participating to the collision ((Apart)) calculated using the Glauber model. 
To account for the bias due to our large centrality bins, (A" part ) has been weighted by the 
number of binary nucleon-nucleon collisions (AT coll ), which is expected to be the scaling 
variable of the J/if> production cross-section in A- A, in absence of nuclear matter effects. 
This correction is small, except for the most peripheral bin where (N paxt ) = 46 while 
the weighted value is 70. These results show no significant dependence on centrality, 
and the integrated R° AA 0% = °- 49 ± 0.03(stat.) ± 0.11(syst.). The comparison with the 
PHENIX measurements at a/snn = 200 GeV [3j H] shows that the inclusive J /if) -Raa 
at 2.76 TeV in the rapidity domain 2.5 < y < 4 are clearly above those measured at 
200 GeV in 1.2 < \y\ < 2.2, while they are closer to the midrapidity values at 200 GeV 
(except in the most central collisions). The contribution from the B feed down to the 
J /if) production in our rapidity and domain has been measured to be ~ 10% in p-p 
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.ALICE (Pb-Pb^s^ = 2.76 TeV), 2.5<y<4, p t >0 (preliminary) 
o PHENIX (Au-Au\/s N ^= 0.2 TeV), 1 .2<|y|<2.2, p r >0 (arXiv:11 03.6269) 
, PHENIX (Au-Au\/s N ^= 0.2 TeV), |y|<0.35, P T >0 (nucl-ex/0611020) 
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Figure 2. Left: J /if) Raa as a function of (N palt ) compared with PHENIX results 
in Au-Au collisions at ^/snn = 200 GeV. Right: J ftp Rcp as a function of centrality 
compared with ATLAS results. Error bars represent the statistical uncertainties, open 
boxes represent the centrality-dependent systematic uncertainties while the centrality 
independent uncertainties are shown by filled boxes. 



collisions at a/snn = 7 TeV [IT] . Therefore, the difference between the prompt J/ip 
Raa and our inclusive measurement is expected to be ~ 11% if the b production scales 
with iV co ii or smaller if it is suppressed by nuclear effects (shadowing, ...). Finally, the 
comparison of our J/ip Rqp results to the ATLAS measurements in the same centrality 
classes [12] (Fig. [fright) indicates that the J/if) mesons measured at forward rapidity 
down to = are less suppressed than the high-p T J/ip mesons at midrapidity (80% 
of the J /if) particles measured by ATLAS have a p<r larger than 6.5 GeV/c). 

In summary, these results show a significant suppression of the inclusive J /if) 
production in Pb-Pb collisions at a/snn = 2.76 TeV. The comparison with PHENIX 
results suggests that re-generation mechanisms could play a role. In order to provide 
tight constraints to suppression/regeneration models, a better knowledge of CNM effects 
is required, which can be obtained with a measurement of J /if) production in p-A 
collisions at the LHC. 
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